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Abstract. To test the competing hypotheses of poly¬ 
phyly and monophyly of “sclerosponges," sequences 
from the 5' end of 28S ribosomal RNA were obtained 
for Astrosclera willeyana, Acantlwcbaetetes wellsi, and 
six other demosponge species. Phylogenetic relationships 
deduced from parsimony and neighbor-joining analyses 
suggest that these sclerosponges belong to two different 
orders of Demospongiae: Astrosclera willeyana, being 
closely related to the Agelasidae, belongs to the Agela- 
sida, Acantlwcbaetetes wellsi, being closely related to the 
Spirastrellidae, belongs to the Hadromerida. These results 
contradict the hypothesis that sclerosponges are mono- 
phyletic and imply that a massive calcareous skeleton has 
evolved independently in several lineages of sponges. 

Introduction 

Recent sponges generally have a skeleton made of spic¬ 
ules that are either siliceous (classes Demospongiae and 
Hexactinellida) or calcareous (class Calcarea). However, 
16 living species build an unusual solid calcareous skele¬ 
ton, which bears a striking similarity to that of various 
Cnidaria, in addition to this spicular skeleton. These 
“coralline sponges” are believed to be the survivors of 
the stromatoporoids, sphinctozoans, and chaetetids, im¬ 
portant ancient reef builders that were highly diversified 
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during the Paleozoic and Mesozoic eras and that were 
long thought to be extinct (Hartman and Goreau, 1970; 
Vacelet, 1983; Wood, 1990). 

Since the discovery of these living coralline sponges, 
they have been classified according to three systems, each 
reflecting a different belief in the number of times that 
sponges have invented a massive calcareous skeleton. In 
the first, all of the coralline sponges are included in the 
class Ischyrospongiae (Termier and Tennier, 1973). In 
the second, the massive calcareous skeleton is believed to 
have evolved at least twice, once among coralline sponges 
with similarities to the Calcarea and once among coralline 
sponges that more closely resemble the Demospongiae; 
the latter group is assigned to the class Sclerospongiae 
{Hartman and Goreau, 1970). This second interpretation 
has been the most widely used, appearing in many recent 
treatises on zoology (Parker, 1982; Riedl, 1983) and pale¬ 
ontology (Rigby and Stearn, 1983). A third system (Va¬ 
celet, 1979, 1985) reflects the assertion that the massive 
calcareous skeleton is more plastic and has evolved in 
several different lineages within the Demospongiae and 
Calcarea. Under this system, living and, where possible, 
fossil coralline sponges are classified within the various 
taxa of Demospongiae and Calcarea with which they 
share derived characters. 

Three coralline sponges are included in this study: 
Acantlwcbaetetes wellsi Hartman and Goreau, 1975; 
Astrosclera willeyana Lister, 1900; and Petrobiona mas- 
siliana Vacelet and Levi, 1958. Acantlwcbaetetes wellsi 
and Astrosclera willeyana are of special interest, because 
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they are considered to be living representatives of chaete- 
tids and stromatoporoids, two groups of great importance 
in the fossil record. The affinities of these groups were 
previously uncertain, but they were most often classified 
in the Cnidaria (Lecompte, 1956; Fischer, 1970). The 
separation of the two groups implies an independent deri¬ 
vation of the massive calcareous skeleton. The spicular 
and cytological characters of both species strongly resem¬ 
ble those found in well-defined families of non-calcified 
demosponges. The choanocytes of Acanthochaetetes 
wellsi possess a periflagellar sleeve; a central cell at the 
apopyle of the choanocyte chambers, as in the Hadromer- 
ida; and a spicule complement similar to that of the family 
Spirastrellidae in the order Hadromerida {Hartman and 
Goreau, 1975; Vacelet and Garrone, 1985; Reitner and 
Engeser, 1987; Boury-Esnault et al., 1990). Astrosclera 
willeyana has small choanocyte chambers, flattened cho¬ 
anocytes, verticillate acanthostyles, and chemical affini¬ 
ties with the order Agelasida ( Hartman and Goreau, 1970; 
Vacelet, 1981; Boury-Esnault et al., 1990; Williams and 
Faulkner, 1996). Petrobiona massiliana has morphologi¬ 
cal affinities with the class Calcarea. 

In this work we generate a new, independent data set 
based on DNA sequences, use it to construct phylogenies, 
and compare these with morphological ones. Our main 
objective is then to determine which of the three hypothe¬ 
ses are consistent with the molecular phylogeny. 

Materials and Methods 

Material: selection and preservation 

The species analysed and their sites of collection are 
listed in Table 1. Some demosponge species were selected 
as representatives of the various taxa supposedly related 
to Astrosclera and Acanthochaetetes. Other species with 
various levels of distance from the in-group were chosen; 
these include representatives of other demosponge sub¬ 
classes—one Ceractinomorpha species ( Halichondria 
panicea) and two Tetractinellida species ( Cinachyrella 
sp. and Discodennia poly discus) —and of class Calcarea 
(Clathrina cerebrum), with which Petrobiona massiliana 
has affinities. For further convenience, all demosponge 
species that are not Tetractinellida are grouped under the 
collective term “monactinellids.” All specimens were 
either preserved in 70% ethanol or deep-frozen in liquid 
nitrogen and then kept at -80°C, depending on collecting 
conditions. 

DNA processing 

Extraction. The total genomic DNA extraction tech¬ 
nique was modified from the Simple Fool’s Guide to 
PCR (Palumbi et al.. 1991). Less than 0.5 g of tissue 


was crushed in a sterile mortar after total dehydration 
(overnight air-dry at +4°C or speed vac) for the alcohol- 
preserved samples and in liquid nitrogen for the frozen 
samples. The powder was gently mixed for a few minutes 
with 500 p\ of lysis buffer (Palumbi etai, 1991). Spicules 
and cellular remains were then removed by centrifugation 
for 2 min at 13,000 rpm. Digested tissue was purified suc¬ 
cessively in phenol, phenol-chloroform-isoamylalcohol, 
and chloroform-isoamylalcohol extractions. Nucleic acids 
were precipitated with ammonium acetate-isopropanol, 
followed by a 70% ethanol wash. Total DNA was resus¬ 
pended in sterile distilled water and its concentration de¬ 
termined by optic density at 260 nm. 

Polymerase chain reaction. Two overlapping frag¬ 
ments of the ribosomal RNA gene were amplified using 
a universal primer and a sponge-specific primer. Primers 
used were as follows (specificity, orientation, and position 
of primers in the aligned sequences of Figure 1 follow 
each sequence); ITS3 5'-GTCGATGAAGAACGCAGC- 
3', universal, forward, external 5'; Eplb' 5'-GTGGCC- 
GGGAGAGGCAGC-3', part of Demospongiae not Tetrac¬ 
tinellida, forward, 257-274; Ep2 5'-CTYYGACGTGCC- 
TTTCCAGGT-3', Demospongiae, reverse, 303-323; D2 
5'-TCCGTGTTTCAAGACGGG-3', universal, reverse, 
external 3'. 

The fragment **ITS3-Ep2" contains a part of the 5.8S 
rRNA gene, the ITS2, the Cl domain and half of the D1 
domain of the 28S rRNA gene; the “Eplb'-D2" fragment 
contains the other half of the D1 domain, the C2 and the 
D2 domains of the 28S rRNA gene. 

A 50 p\ double-stranded PCR reaction mix contains 
0.3 pg template DNA, 2.5 p\ DMSO, 0.165 mA/ each 
dNTP, 30pmol each probe, 1.5 U Taq DNA polymerase 
(Bioprobe). This reaction mix was overlaid with mineral 
oil and placed in a Trio-thermoblock thermocycler (Bio- 
metra). Cycling conditions are variable for the annealing 
temperature (Ta): respectively 60°C and 63°C for ITS3- 
Ep2 and Eplb'-D2 primer pair. The first cycle is 4 min 
at 94°C, 2 min at Ta, and 2 min at 72°C; this is followed 
by 30 cycles each consisting of 1 min at 94°C, 1 min at 
Ta, and 1 min at 72°C; the reaction is finished by 4 min 
at 72°C. 

After visualization of 5 p] of the reaction on a 1.5% 
agarose gel, the remaining 45 pi of PCR product was 
purified by precipitation with ammonium acetate-isopro¬ 
panol, followed by a 70% ethanol wash. The pellet was 
then resuspended in 6 p\ of sterile distilled water. 

The approximate concentration was evaluated visually 
by electrophoresis of 1 pi of the purified PCR product in 
a 1.5% agarose gel, and comparison to 1.5 pi of the DNA 
molecular weight marker VI (Boehringer Mannheim). 

Cloning and sequencing. Each PCR fragment was 
cloned into PCR-Script SK( + ) cloning vector (PCR- 
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Table I 


Sponge species sequenced for analysis of phylogenetic relationships among sclerosponges 


Classification 

Species 

Collection locality 

DEMOSPONGIAE 

Tetractineliida 

Tetillidae 

Cinachyrella sp.* 

New Caledonia 

Theonellidae 

Discodermia polydiscus Bocage, 1870* 

Mediterranean sea, 3PP cave. La Ciotat 

“monactinellids” 

Axinellidae 

Axmella damicornis (Esper, 1794)* 

Mediterranean sea. La Ciotat 

Agelasidae 

Agelas oroides (Schmidt, 1864)** 

Mediterranean sea. La Ciotat 

Astroscleridae 

Astroscleral willeyana Lister, 1900 

New Caledonia 1992 


Astroscleral willeyana Lister, 1900 

New Caledonia 1994 

Clionidae 

Cliona viridis (Schmidt, 1862)* 

Mediterranean sea. La Vesse 

Spirastrellidae 

Spirastrella cf coccinea (Duchassaing & Michelotti, 1874) 

Panama, Atlantic coast San Bias Island 

Acanthochaetetidae 

Acanthochaetetes wellsi Hartman & Goreau, 1975 

New Caledonia 1992 

Halichondriidae 

Halichondria panicea (Pallas, 1766)* 

South West Channel, Aber Wrac’h 

CALCAREA 

Clathrinidae 

Clathrina cerebrum (Haeckel, 1872)** 

Mediterranean sea. La Vesse 

Petrobionidae 

Petrobiona massiliana Vacelet & Levi. 1958** 

Mediterranean sea. Anse des Cuivres 


* Sequences from Chombard et al. (In press). 
** Sequences from Lafay et al. (1992). 


Script SK(+) cloning kit. Stratagene) and sequenced with 
the T7 Sequencing kit (Pharmacia Biotech) using [33P]- 
dATP and adding DMSO in ihe annealing reaction. The 
internal probe C2' is used to obtain the middle of the 
“Eplb'-D2” fragments, in addition to the vector probes 
Ks and T3 (C2' 5'-GAAAAGAACTTTGRARAG- 
AGAGT-3', universal specificity, forward orientation, po¬ 
sition 483-505 on the aligned sequences of Figure 1). 

Each PCR product was sequenced from a minimum of 
two clones; when contradictions in the sequences of sev¬ 
eral clones could not be resolved, the corresponding posi¬ 
tions were coded according to the UP1AC code. The two 
strands were sequenced for the main part of the sequence 
length, with special attention to the D2 domain where 
strong secondary structures of the molecule cause com¬ 
pressions in the sequence migration. From the (wo over¬ 
lapping PCR products, the final sequence was 1104 bp to 
1197 bp in length, depending on the species. This frag¬ 
ment corresponds to the 3' extremity of the 5.8S rRNA 
(about 108 bp), the Internal Transcribe Spacer ITS2 (be¬ 
tween 167 bp and 224 bp), and the four first domains of 
the 5' extremity of the 28S rRNA: Cl, Dl, C2, and D2 
(between 816 and 866 bp). 

Sequence management and alignment 

The MUST package (Philippe, 1993) was used to man¬ 
age sequences, including registration (with ENTRYSEQ 
program), alignment (with ED), construction of distance 


matrices (with NET or from NJ trees), distance calcula¬ 
tions and construction of trees with the neighbor joining 
algorithm (with NJ), matrix comparison (with COMP- 
MAT), and calculation of bootstrap proportions from 
neighbor joining trees (with NJBOOT). Wherever likely 
secondary structures were detected, sequences were 
aligned according to supposed conservation of helices. 

PAUP, version 3.1.1 (Swofford, 1991), was also used 
for construction of trees and calculation of bootstrap pro¬ 
portions, discussed below. In bootstrap calculations, non¬ 
majority nodes were compared in order to explore the 
robustness of alternative topologies. 

The final alignment presented in Figure 1 was obtained 
by eye using the editor of MUST (ED). The 1TS2 (not 
presented in Fig. 1) and part of the 5' extremity of the 
D2 domain (corresponding to positions 575-640, Fig. 1) 
are very divergent and cannot be aligned in all our sam¬ 
ples, thus these regions were not used in the sequence 
analysis. 

Results 

Because previously published sequences of 28S rRNA 
(Lafay et al., 1992) are shorter than ours, two successive 
analyses were made. The first grouped all species and 
corresponds to the length published by Lafay et al. (Table 
I); in the second, Clathrina, Petrobiona, and Agelas were 
removed so that we could use our total alignable length. 

The first analysis included 12 species and 374 bp of 
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Astrosclsral willeyana 
Aatroaclera2 willeyana 
Axinella damicomie 
Acanthochaetetes wellsi 
Spiraatrella cf. coccinea 
Cliona viridia 
Halichondria panicea 
Discodermia polydiecua 
Cinachyrella ap. 


1 10 20 30 40 50 60 70 80 90 

AAACTGCGATACGTAGTGTGAATTGCAGAATTCAGTGAATCATCGAGTCTTTGAACGCAAATGGCGCTCCCGGTCAAGCCGGGAGCACGT 


c- 


--C- 
-T-- 


G-c- 

G- 


c- 


-G- 


Y- 


-T---. 

-T-TTCA- 

-T-TTCA- 

-T-TTCA- 

--T-T-CT-- 

-CGT- 

TCGT-GG--- 


-A- 

-—C*-S-- 

—c*- 


-s- 


Astroscleral willeyana 
Aatroaclera2 willeyana 
Agelaa oroidea 
Axinella damicomia 
Acanthochaatatas wellai 
Spiraatralla cf. coccinaa 
Cliona viridia 
Halichondria panicea 
Diecodermia polydiecua 
Cinachyrella ap. 

Clathrina cerebrum 
Petrobiona maaailiana 


5.8S][28S-C1 

100 110 120 130 140 150 160 170 180 

CTGT CT gagcgtccttttttggacctcagatcaggc gaggct acctgc tgaactt aagcatatcaat aagcagaggaaaagaaact aac a 

-vf- 

NN-A- 

-TG-ATC-T- 

-A-TTG-C—A-T-ATC-C-T-G- 

-A-T-GA--C--A-T-Y-- ATC-C-T-G- 

-A-T-GA—C—A-T- ATC-C-T- G- 

-A--TTGC--C--A-A-C-T-G- 

-G--GC-G-C-G-C-C 

-CA-G-C-G-C-C 

NNN -G-T-A--A-C-T-T-G- 

NNT -G-T-A-C-T-»**-G- 


Aatroacleral willeyana 
Aatroaclera2 willeyana 
Agalaa oroidea 
Axinella damicomia 
Acanthochaatetea wellai 
Spiraatrella cf. coccinea 
Cliona viridia 
Halichondria panicea 
Diacodarmia polydiacua 
Cinachyrella ap. 

Clathrina cerebrum 
Petrobiona maaailiana 


Cl][Dl 

190 200 210 220 230 240 250 260 270 

AGGATTC CCCCAGT AAC GGCGAGC GAAGTGGGAAAAGCTCGAGC C TGAAATCT CTGGCAGTTGATGCCAGCGAATTGTGGC CGGGAGAGG 


■T- 


-T-- 

--T- 


•T- 


-C-G-A-T- 

- C -CG—--C-G--T-A-A- 

- c - c .G--C-A- 

-c -C-C-G--C-A- 

- G - c -CAC--T-G—G- 

-C-C-AG-C-TGTCGATC-C- 

-C-TG-AGTC-GTGGAGC- 

-T-C-C-ATT-A-GC--TGC-CAT-GTCA-TTTAA-A- 

■ -T-C-G-A-GA-GC-CGT-TTCA-N-TT-A-A- 


Aetroaclaral willeyana 
Aetroaclera2 willeyana 
Agalaa oroidea 
Axinalla damicomia 
Acanthocheetetea wellai 
Spiraatrella cf. coccinea 
Cliona viridia 
Halichondria panicea 
Diacodermia polydiacua 
Cinachyrella ap. 

Clathrina cerebrum 
Petrobiona maaailiana 


280 290 300 310 320 330 340 350 360 

CAGCTGGACCCTGGCT AGCGCTGTCGAAGTTGACCTGGAAAGGCACGTCAGAGAGGGTGACAGCCCCGTGTGCGACACTGC * CGGGCAGG 

-G-*- 

-CG-C--*-* *- 

-CG-T-A-CA-*T—A- 

-G-TG-G-A--GTT-A--G-G-CTT-C-A*-A--- 

-G--G-TG-G-A--GTT-A-G-C-A*-A- 

-G—GTTG-CTA—GT-G-G-C-C-A*-A 

--A--A-G-A-CGA-A-G-C-G-AC-T-G-C-T*--AC- 

-G-G--C-AGACGCGA-T-C-G-T-C G* CGGAGT - -CCG-A- 

-G-G--C--TGGA-GCGC-T-G-T-G-T-GGCGGAGGT-CCG-AA 

TGTT - TCC - -GG - ATGTC-G-CT-AG-NG-T-CA-T-N--G-A-C--T-G-T-G-TGGANN-TT- 

TGCT-TC-GG-A-G-CG-A—G-CT-AG-NG-T-CA-T-G-AT-C--T-GTTGGT-GGTC—TC- 


Astroaclaral willeyana 
Aatroaclara2 willeyana 
Agalaa oroidea 
Axinella damicomia 
Acanthochaetetea wellai 
Spiraatrella cf. coccinea 
Cliona viridia 
Halichondria panicea 
Diacodermia polydiacua 
Cinachyrella ap. 

Clathrina cerebrum 
Petrohiona maaailiana 


Dl][C2 

370 380 390 400 410 420 430 440 450 

GCCACCACTGTCTTCAGAGAGTCGGGTTGTTTGGGAATGCAGCCCAAAGTGGGTGGTAAACTCCATCTAAAGCTAAATATGGGCACGAGA 


-N-. 

-T-G-T- 

-G-T-G— 

-G-T—C-G-- 

-G-TT-C-T-- 

-GT--T--C- 

---G-T—C-C-G-- 

--MG-T--C-C-- 

-TATGATGCACT-GA- 

-NATGATGCACT-GA- 


- A __ T -TG- 

-T--T-- 

•N-NT- -A- - 


■N- 

- G - CA - T . 

-CC- 

- CCT - 

-cc- 

-R-CT- 

- C -AC- 

-c-— C-GC- 

-T-C- 

-T-C- 


Figure 1 . Aligned sequences. Only nucleotides that differ from those of Astroscleral are indicated 
(identities are noted by hyphens and deletions by stars). Boundaries between 5.8S gene and 28S gene are 
indicated over the sequences, as boundaries between domains of the 28S gene. Crosses over sequences 
indicate the nonalignable part of the D2 domain, which is not used for phylogenetic analysis. 


sequence, 145 of which are variable and 106 informative 
for parsimony. As shown in Figure 1, these sequences 
include the Cl, Dl, and part of the C2 domains of the 28S 


rRNA gene. Saturation was tested using COMP-MAT 
of MUST. Global saturation is not detected, observed 
distances and number of steps inferred by PAUP between 
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Aetroscleral willayane 
Astroeclara2 willayane 
Agalae oroides 
Axinella damicomis 
Acanthocheetatas wellei 
Spirastrella cf. coccinaa 
Cliona viridis 
Hslichondria panicea 
Discodermia polydiscus 
Cinachyrella sp. 

Clathrina cerebrum 
Petrohiona massiliana 


C2][D2 

460 470 480 490 500 510 520 530 540 

CCGATAGCAAACAAGTACCGTGAGGGAAAGGTGAAAAGTACTTTGAAAAGAGAGTCAAAAAGACCGCGAAACCGTTAGGAGGGAAGCGAA 


-G- 

-G- 

-G- 

-G- 

-G- 


-C-N 


-A- 

-A- 

-A- 


-G- 

G-A... 

-G---A 


-C- 

-c- 


-R- 


■G-T-G-G- 

-G-T-G- 


-A- 

-A- 

A- 

-A- 

-A--C- 
-A—C- 


550 560 570 500 590 600 610 620 630 

Astrosclerel willayane TGCAGCCAAAGTGGTTCTCGTTCAGGCTCAGGAG**TTGTTGGCGTGCAGTGCTGGGATGCCAGACGCCGTAGGGTGCTGCAACTCGGAT 

Aatroeclera2 willeyana -**- 

Axinella damicornis ---C--C-A-* *-C-C---C--TTG-.A-G-CA- 

Acenthocheetetes wellei -TG-C-TCGTA-T.G-* *AC--C--A-GA-*-TC--TG-T-GG-CA-C-*C--C-TCG*C-* 

Spirastrella cf. coccinea-TG-C--CGCA-T-G-**AC--C--GAGAA*-C-TC--C*--GTG-CA-—C-'C-TC-TCGG-G* 

Cliona viridis -TG-C--CGTA-T-G-**AC--C--A-G******TC-ATC--GGTC-A--A--T*C-A*--CCGGA-CG* 

Halichondrie panicea -TG--T-A-C-GCC*--A-T-CGAC-ATCA-TGA-C--TCC--A-A.ATA-TT-CCTC- 

Discodermia polydiscus -GG-C C* -CGCG--GA-T-G-*GCC-CC—GCGT-CCGCTC-**TC*GTC---A--C-C-C*G-G-CGTCGG-*-GG 

Cinachyrella ap. -G--C C*--TCG--G--T C-G-*G-C-C*-CGAGT-CTG-TC-A-CC-GGC-T-A--CAT-C-G-G-CGTCTG-TCGG 


Aatroacleral willeyena 
Astrosclara2 willayene 
Axinella damicornis 
Acanthochaetatas wellei 
Spiraatralle cf. coccinea 
Cliona viridie 
Halichondria panicee 
Diecodermie polydiscus 
Cinachyrella sp. 


640 650 660 

ACGGCTGT** CGACTGCTTTGCATTCCTGACGAGAG 


670 680 690 700 710 720 

**CCGGCCAACGGCAGTTA*CCCCTGGCTCAAGAGGGTTGTTGGGAAGGTAGC 


—A- -A-C* *-C- 

-* TCT-G*-CAG-C- 

-TCCCG*-AA--C- 

G-T*»***-AA--C- 

-**-* * -A-C--C- 

TG-CGGCGTCG-G--C-CC- 

CG-AGGTGTC-*G--C-CC- 


-T-G-* * *T- 

-G-TG-G--- 

-C-TG-G ***- 

- C -TT-G ***- 

-C—T — GGC-*— * * * -T-R-- 

-C--TC-G--GC--GGG- 

- c -C-G--GTC-GGG- 


-A-*-T--A-GG--T--CC-G- 

-G-A-G*T-C--GA-TCG--C-*-GC- 

-G*-A--C--GA-TCG--C-*-GC- 

-A-G*TTT-C-C--GA-TCG--C-*-GC- 

-AA-GC-TTG-G-CA-CC--C---T-GC- 

-GT--G-CGGATT--G*-GCG-CACC--C C- 

-G-A-G--GGTGT--GC-GCG-ACG--C T-G-T 


Astroeclaral willeyana 
Aetrosclera2 willeyana 
Axinelle damicomis 
Acanthochaetatas walls! 
Spirastrella cf. coccinaa 
Cliona viridia 
Halichondria penicaa 
Discodermia polydiscus 
Cinachyrella ep. 


730 

TTCTCGGTT* 

-T-* 

-CTG-AC-* * 
-CTG-AC-** 
CCTG-AC-** 

CA-CT*** 

GCAG-C--**! 


740 750 760 770 780 790 800 

* TACCGGGAAGAACTTACAGCCGGCAACCTGGCA * GTCTGGGAGTGACTGAGGAGTGCTGTGACT * 

* -a-G-ATTGG-*-*-G-C-R—*’ 

* *CGT-CAG G-G--G-T-CGT-C—TG*-GAC-A-TC-A-*’ 

* *CGT-CAG-G-T-CGT-C-G*-C-GAC-A-TC-A- * 

* *CGT-CAGG--G-G--T-TTCGT-C — TG*-C-A--C-A- * 

* * * -G-A-TG-TG-T-GGT-A-*-CAG-TAA*** ****** 

G**-G-A-GC-T-GGTGC--AGC-CTC--CGTCCG-C-A--G-G-A* 


810 

TTTCA 


GCAG-C-GCCCTCTCGGGT-G-A-GC--G- 


-T-CGTTC-G-GC--GTGC--C-A--G-G--TTTG- 


Astroscleral willeyana 
Astroeclera2 willeyana 
Axinella damicoraie 
Acanthocheetatas wellei 
Spirastrella cf. coccinea 
Cliona viridis 
Halichondrie panicea 
Discodermia polydiecus 
Cinachyrella sp. 


820 830 

CACCCGCAGTACAGGCTCCCT* 


840 850 860 870 880 

*AGGGGGCC *GGGTCCCTTCTGTCTGTGTGGGCAACGCCGCAGGG 


890 900 

ACTGCAT* GCAGTGTCTG 


-A-G-******-*__*-G-C-T--TT-T-*-G-*-G- 

-G-GT-CG-CTG-CTC* * *GT-C-*-G-C-C--T-CC-G-TG-*-T-*- 

-G-GT-CG-C-T**‘*******C-*-G-C-C--T CC-G-TG *-T-T-C- 

-G-GT-CG-c--T********C-*-C-A-C-C--T TCC-G-TG *-G-*-C- 

* *-TG-CGT-c--T********T-TTT-C-G-G-C--T-C-TT-*-G-*- 

-C-GT-CG-TC * -GCCT* *CGGC- *-*-G-G--CTC-C--TC--CG-C-CGAC-GC-TC-C*-C- 

-C—G-CG-CT-GCTTTAGGGC-*—AG-G-C—TCA-C-TC-GC-GA—G-TC*-C-C- 


Astrosclerel willayane 
Astrosclera2 willayane 
Axinalla damicomis 
Acenthochaatetee wellsi 
Spirastrella cf. coccinea 
Cliona viridie 
Halichondria panicee 
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Figure 1. (Continued) 


all the pairs of species in the data set being linearly corre¬ 
lated (CC 0.98, Fig. 2). Three groups of dots are clearly 
detectable in this saturation analysis. They correspond to 


decreasing distances between (1) Calcarea-Demospon- 
giae, (2) Tetractinellida-monactinellids, and (3) monacti- 
nellids-monactinellids. The exhaustive search algorithm 
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Pairwise number of inferred substitutions 
in the most parsimonious tree 

Figure 2. Global saturation curve for 12 species and short-length 
aligned sequences (374 bp). CC 0.98. White circles are distances be¬ 
tween pairs monactinellids-monactinellids. Dark circles are distances 
between pairs Tetractinellida-monactinellids. Dark squares are distances 
between pairs Calcarea-Demospongiae. White squares are distances 
within Tetraclinellida and within Calcarea. 


of PAUP provided one single shortest tree with 268 steps, 
a consistency index (Cl) of 0.795, a retention index (RI) 
of 0.767, and a G1 of -1.50. The tree was rooted using the 
out-group method on both species of Calcarea ( Clathrina 
cerebrum and Petrobiona niassiliana). The resulting sin¬ 
gle topology is presented in Figure 3. The Branch and 
Bound search option was used to provide a bootstrap with 
1000 replicates in PAUP. The majority-rule consensus 



Astrosclera 1 willeyana 
Astrosclera2 willeyana 
Age las oroides 
Axinella damicomis 
Acanthochaeteles wellsi 
Spirastrella cf. coccinea 
Cliona viridis 
Halichondria panicea 
Discodermia polydiscus 
Cinachyrella sp. 
Clathrina cerebrum 
Petrobiona massiliana 


Figure 3. Phylogram obtained with PAUP by exhaustive analysis on 
short-length aligned sequences (374 bp) for 12 species using ACCTRAN 
optimization option. Tree length = 268, Cl = 0.795. RI = 0.767, and 
G1 = - 1.50. Bootstrap proportions (1000 replicates using Branch and 
Bound) are shown above internal branches. 



Astrosclera 1 willeyana 
Astrosclera2 willeyana 
Agelas oroides 
Axinella damicomis 
Halichondria panicea 
Acanthochaeteles wellsi 
Spirastrella cf. coccinea 
Cliona viridis 
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. Clathrina cerebrum 
. Petrobiona massiliana 


Figure 4. Phylogram obtained with MUST by neighbor-joining 
analysis on short-length aligned sequences (374 bp) for 12 species. Boot¬ 
strap proportions (1000 replicates by NJ analysis) are shown above 
internal branches. 


tree exhibits the same topology as the shortest tree found 
by exhaustive search (bootstrap proportions IBP] are re¬ 
ported on Fig. 3). The neighbor joining analysis (NJ and 
NJBOOT in MUST) provided a topology that differs in 
the location of Halichondria and of Acanthochaetetes, 
and in having slightly better bootstrap proportions 
(Fig. 4). 

This first analysis indicates that the Tetractinellida are 
monophyletic, a conclusion supported by a 100# BP 
(Chombard et al. , in press). This group constitutes the 
sister group of the other Demospongiae called here “mo- 
nactinellids." This last group is supported by a 96# BP 
in distance analysis and an 84# BP in parsimony analysis 
(Figs. 3-4). All the alternative topologies found by parsi¬ 
mony have less than 5# BP. implying that monactinellids 
are the monophyletic sister group to the Tetractinellida. 
For the second analysis of full-length sequences, we are 
thus able to take the Tetractinellida as an out-group re¬ 
lated to the monactinellids. In monactinellids, “sclero- 
sponges" are polyphyletic. Acanthochaetetes is included 
in a hadromerid clade, in which the monophyly of (Acan¬ 
thochaetetes, Spirastrella, Cliona) is supported by respec¬ 
tively 100# BP in distance and 99# BP in parsimony 
analysis (Figs. 3-4). Relationships within this clade are 
not strongly supported by this first analysis. Astrosclera 
(two individuals) is included in an axinellid clade, in 
which the monophyly of (Astrosclera 1 , Astrosclera2, 
Agelas, Axinella ) is supported by 96# BP and 83# BP 
in distance and parsimony analysis respectively. Unlike 
the hadromerid clade, the axinellid clade has relationships 
that are well supported—in particular the monophyly of 
(Agelas, Astrosclera1, Astrosclcra2 ), which is supported 
by 100# BP and 99# BP in distance and parsimony anal¬ 
ysis respectively. 

The second analysis was made for 9 species and 914 bp 
of sequence, 388 of which are variable and 244 informa- 
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Pairwise number of inferred substitutions 
in the most parsimonious tree 

Figure 5. Global saturation curve for 9 species and full-length 
aligned sequences (914 bp). CC 0.97. White circles are distances be¬ 
tween pairs monactinellids-monactinellids. Dark circles are distances 
between pairs Tetractinellida-monactinellids. While squares are dis¬ 
tances within Tetractinellida. 


live for parsimony. No global saturation is evident (CC 
0.97, Fig. 5). The exhaustive search algorithm of PAUP 
provided one single shortest tree with 640 steps. Cl = 
0.839, R1 = 0.783, and G1 = -1.01. The iree was rooted 
using ihe out-group method on the tetractinellids (Cina- 
chyrella sp. and Discodennia polydiscus). The resulting 
single topology is presenied in Figure 6. The Branch and 
Bound search option was used to provide a bootstrap with 
1000 replicates in PAUP. The majority-rule consensus 
tree exhibits ihe same topology as the shortest tree found 
by exhauslive search (BP reported on Fig. 6). The neigh- 
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Figure 6. Phylograni obtained with PAUP by exhaustive analysis on 
full-length aligned sequences (914 bp) for 9 species using ACCTRAN 
optimization option. Tree length = 640. Cl = 0.839, R1 = 0.783, and 
G1 = —1.01. Bootstrap proportions (1000 replicates using Branch and 
Bound) are shown above internal branches. 
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Figure 7. Phylograni obtained with MUST by neighbor-joining 
analysis on full-length aligned sequences (914 bp) for 9 species. Bool- 
strap proportions (1000 replicates by NJ analysis) are shown above 
internal branches. 


bor-joining analysis (NJ and NJBOOT in MUST) pro¬ 
vided the same topology and similar bootstrap proportions 
(Fig. 7). This second analysis confirms the first one: ihe 
sclerosponges Astrosclera and Acanthochaetetes belong 
to two different clades, a hadromerid clade and an axinel- 
lid clade. The hadromerid clade (Cliona, Spirastrella, 
Acanthochaetetes) is supported by 100% BP in both dis¬ 
tance and parsimony analysis, and the internal topology is 
also supported by 92% BP and 88% BP for (Spirastrella, 
Acanthochaetetes) in dislance and parsimony analysis. 
The axinellid clade ( Axinella, Astrosclera], Astrosclera2) 
is supported by 100% BP. The monophyly of the two 
Astrosclera individuals is supported by 100% BP; ihe in¬ 
dividuals came from the same area of New Caledonia and 
do nol represent Ihe two populations, differing by the 
presence or absence of spicules, thal occur respectively 
in the Indian Ocean and the Central Pacific (Vacelet, 
1981; Ayling, 1982). 


Discussion 

The Ischyrospongiae (Termier and Termier, 1973) hy- 
polhesis is falsified by ihe first analysis. The “coralline” 
sponge Petrobiona inassiliana clearly belongs to the class 
Calcarea, whereas the two other calcified sponges, Astro¬ 
sclera and Acanthochaetetes, are undoubtedly part of the 
Demospongiae. The class Ischyrospongiae is thus poly- 
phyletic, as concluded previously from morphology, and 
should be abandoned. 

Both current analyses demonslrate ihe polyphyly of (he 
class Sclerospongiae, and it too should be abandoned in 
classification schemes. Furthermore, the two sclero¬ 
sponges belong to different monophyletic clades, an axi¬ 
nellid one ( Axinella , Agelas, Astrosclera) and a hadrom¬ 
erid one ( Cliona, Spirastrella. Acanthochaetetes). Both 
clades are strongly supported in the iwo analyses. They 
are in complete agreement with the affinities indicated 
by spicule morphology and by cytology (Hartman and 
Goreau, 1970, 1975; Vacelet, 1981; Vacelet and Garrone, 
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1985; Reitner and Engeser, 1987; Boury-Esnault et al .. 

1990) . These results support the interpretation that the 
capacity to secrete a massive skeleton of calcium carbon¬ 
ate has developed several times during the course of the 
evolution of the Porifera (Vacelet, 1979, 1983, 1985; 
Wood et al., 1989). Accordingly, these “coralline" 
sponges have to be classified in the Demospongiae; Acan- 
thochaetetes wellsi in the order Hadromerida; and Astro- 
sclera willeyana in the order Agelasida, which is consid¬ 
ered by most recent authors as distinct from the order 
Axinellida, although closely related to it. The creation 
of a special order- -the Tabulospongida—based on the 
presence of a calcareous skeleton (Hartman and Goreau, 
1975) in Acanthoehaetetes wellsi and its fossil relatives 
has no strong justification according to the present results. 

At a lower taxonomic level, the classification of these 
sponges as belonging either within existing families to 
which they are closely related or in distinct families is 
still subjective. Pending analyses of other related sponges, 
the decision depends upon individual judgments about 
the size of the morphological gap needed to separate taxa 
and about the importance of the calcareous skeleton as a 
taxonomic character. In the case of Acanthoehaetetes, we 
propose to classify the genus in the family Spirastrellidae 
Ridley and Dendy, 1886, in view of the spicular and 
cytological resemblances (periflagellar sleeve, central 
cell) and the low genetic distance between Acantho- 
chaetetes and Spirastrella that is indicated by the present 
work. This hypothesis, which was already proposed by 
Reitner (1991), avoids the use of the family Acantho- 
chaetetidae, which would be monogeneric at least in the 
Recent fauna. (We reject, however, on the grounds of 
morphology, Reitner's merging of the genus Acantho- 
chaetetes with Spirastrella.) In the case of Astrosclera, 
we prefer to maintain the two families Astroscleridae 
Lister, 1900 (with five genera in the Recent, if merged 
with Ceratoporellidae), and Agelasidae Vcrril, 1907 (with 
one large genus). The genetic distance between Astro¬ 
sclera and Agelas, as estimated by our sequences, is ad¬ 
mittedly as low as for Acanthoehaetetes and Spirastrella. 
However, the Agelasidae and the Astroscleridae differ by 
an important reproductive character: Agelas is oviparous 
(Liaci and Sciscioli, 1975; Reiswig, 1976), whereas 
Astrosclera is viviparous (Lister, 1900). Furthermore, the 
structure of the spongin fibers of Agelas (De Vos et al., 

1991) , which are unique among the Demospongiae, is an 
important difference between Agelas and Astrosclera. 
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